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Abstract 

Motivated by the paper |SWj by B.T. Sutcliffe and R.G. Woolley, we present the 
main ideas used by mathematicians to show the accuracy of the Born-Oppenheimer 
approximation for molecules. Based on mathematical works on this approximation 
for molecular bound states, in scattering theory, in resonance theory, and for short 
time evolution, we give an overview of rigourous results obtained up to now. We also 
point out the main difficulties mathematicians are trying to overcome and speculate 
on further developments. The mathematical approach seems not to fit exactly to 
the common use of the approximation in Chemistry. We contribute in this way to 
the discussion on the Born-Oppenheimer approximation initiated in |SW ] . 
The paper neither contains mathematical statements nor proofs. Instead we try to 
make accessible mathematically rigourous results on the subject to researchers in 
Quantum Chemistry or Physics. 



1 Introduction. 



In the paper [SWJ . the authors made a remarkable effort to understand the mathematical 
literature on the Born-Oppenheimer approximation. It was certainly not a easy task for 
them to extract relevant information for Chemistry from papers, which often use elaborate 
mathematical tools and provide more or less abstract results. For instance, they comment 
on the paper |KMSW] , that makes use of semiclassical pseudodifferential calculus and of an 
important, but rather complicated trick (due to Hunziker in [Huj ) to control the Coulomb 
singularities appearing in the potential energy of the molecule. They also pointed out to 
their colleagues in Chemistry some misunderstandings and too crude simplications in the 
traditional treatment of the Born-Oppenheimer approximation for molecules. One may 
feel a slightly pessimistic note in the paper |SW] on the possibility for Chemists to use 
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the Born-Oppenheimer approximation in a correct and accurate way and to benefit from 
mathematical works on the subject. Here we shall give a description of the situation from 
the mathematical point of view, which does show the present difficulties and limitations 
of the mathematical approach but still leads to a quite optimistic impression. This is 
probably due to a slightly different interpretation of the works |B0] and (BHj . 
In the last three decades, mathematically rigourous works on the validity of the Born- 
Oppenheimer approximation for molecules has been produced. We sort these works in 
two (uncomplete) lists into alphabetic order. The first one contains articles that, strictly 
speaking, study the Born-Oppenheimer approximation: [Cl ICDS| ICSl IHal' 'Ha21 IHa3[ 

[M iHi5i ihtti iHHi [HjTi nmmmi imni imsi ihrji in i \m \M ijkwi ikmsw. kmwh 

IKMW21 Man Ma2[ MaSi MM iMMl IMSTI [MS2l IPSTI Eil El |ST1 |TW]. In the second 
list, we mention closely related works on semi classical Schrodinger matrixoperators: |DFJl 
WUi [FRl IFLNI El [Hi6l [HJ4l [H |HT]. In the first list, the papers essentially show that a 
reduced Hamiltonian (a Schrodinger operator with matrix or operator valued potential) 
is a good approximation to the true molecular Hamiltonian. In the second list, the works 
obtain mathematical results on the reduced Hamiltonian, that are of physical or chemical 
relevance for molecules. 

In the present paper, we focus on the Born-Oppenheimer approximation, namely the 
possibility to approximate the true molecular Hamiltonian by some effective Hamiltonian 
usually called the adiabatic operator. In Section [2|, we introduce the Hamiltonian of 
the system and proceed to the removal of the centre of mass motion in two ways, one 
adapted to the study of bounds state and the other to scattering theory. In Section [3], we 
present the core of the mathematical form of the Born-Oppenheimer approximation and 
describe the construction of the adiabatic Hamiltonian. In Section IH we explain selected 
mathematical works and comment on the actual difficulties and limitations of the theory. 
Finally, in Section |5l we sum up the main features in the mathematical Born-Oppenheimer 
approximation and argue that further progress towards chemically relevant questions can 
be reasonably achieved. Two figures used at many places in the text are added at the end 
of the paper. 

As pointed out in the abstract, we only present intuitive arguments and statements, that 
do not respect at all the standard rigour in mathematics. But they do have a rigourous 
counterpart in the mathematical literature. 

Acknowledgement: The author is particularly grateful to S. Golenia who performs the 
integration of the figures into the source file of the text. He also thanks B.T. Sutcliffe and 
R.G. WooUey for fruitful discussions. The author is financially supported by the CNRS 
and the ANR "NOSEVOL". 



2 The Hamiltonian. 

We consider a molecule with M nuclei with positive masses mi, m2, ■ ■ ■ , mM respectively, 
with positive charges Zi,Z2,- ■■ ,Zm respectively, and with electrons with mass set 
equal to 1. We set the Planck constant h and the electronic charge e to 1. Denoting by 
zi,Z2, - ■ ■ , zm G M'^ the positions of the nuclei and by zm+i, zm+2, ■ ■ ■ , zm+n £ the 
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positions of the electrons, the Hamiltonian of the molecule is given by: 

Pmoi = K + W, (2.1) 

M ^ ^ M+N 

k=l j=M+l 

jjG{l,-,Af} ' i,ie{Af+l,---,J\/+Ar} ' ' 

M M+N „ 

-EE ■ (2-3) 



Here — A^^. denotes the Laplace operator in the = {zl, zl, zf) variable, that is 

where dt stands for the partial derivative with respect to the variable t. 

It is usual and physically relevant to remove from the Hamiltonian Pmoi the motion of 
the centre of mass of the molecule. This is done by an appropriate change of variables. 
There is no canonical choice for this change of variables. This means in particular that 
one can choose it according to the kind of study one wants to perform. To study bound 
states of the molecule or its time evolution, we shall use the change of variables adopted 
in [KMSWl IMS21 ISW] . To consider diatomic collisions (ion- ion, ion-atom, or atom-atom 
scattering), we shall use another one (those in |KMWH IKMW2j ): see p. 75-82 in |RS3j 
for details on the removal of the centre of mass. 

In the first mentioned situation, we take the nuclear centre of mass (which is close to the 
centre of mass of the molecule), Jacobi coordinates for the nuclei, and atomic coordinates 
for the electrons. Let C : (zi, ■ ■ ■ , zm+n) ^ {R] a^i, ■ ■ ■ , xm-i, Hi, - ■ ■ , Vn) be the change 
of variables defined by 



M ^ M 



m 

k=l k=l 



fori < j < M - 1 , Xj = Zj+i - nikZk 



R is the centre of mass of the nuclei, the xj are the new "nuclear" coordinates, and the 
yj are the new electronic variables. For an appropriate constant C that only depends on 
the masses and on N we define, for any L^ function / of the variables {zi, ■ ■ ■ , zm+n), 

{Uf){R;xi,- ■ ■ ,XM-i,yi,- ■ ■ ^Vn) = C/(C~^(i?;xi, • • ■ ,xm-i,2/i, ■ ■ ■ ,l/7v)) • (2-4) 

The constant C is chosen such that, for all /, / and Uf have the same L^-norm {U is 
unitary), keeping unchanged the physical interpretation of the L^-norm. Looking at the 
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Hamiltonian Pmoi in the new variables {R; Xi, ■ ■ ■ xm-i, Vir ' ' ! Vn) means that we consider 
the operator 

Im 

where H only acts on the variables (xi, ■ ■ ■ ,XM-i,yi, ■ ' ' ^Vn)- Forgetting about the 
kinetic operator of the nuclear centre of mass, we focus on the physically relevant Hamil- 
tonian H. Denoting by Vt the gradient operator in the variable t and setting, 

for 1 < j < M- 1 , fij^ = mjl^ + rrik^ , x = {xi, ■ ■ ■ , Xm-i) , and y = [yi, ■ ■ ■ , i/at) , 
the latter is given by 

H = Ho + TuE, (2.5) 

M-l ^ 

The = ^k,Vy,-Vy^, (2.7) 

^<k<j<N 
N 

Q{x) = -$^A,, + W{x-y) (2.8) 

k=l 

where W{x; y) is just the function W in (12. 3p composed with the inverse change of variables 
C^^. We observe that Q{x) is an operator in the electronic y variables that depends only 
parametrically on the nuclear x variables and does not depend on R. The operator The 
is usually called the Hughes-Eckart term. For each nuclear configuration x, the operator 
Q{x) is refered to as the electronic Hamiltonian in the configuration x (it is called the 
clamped- nuclei Hamiltonian in |SWj ) . Note that the coefficients fij in H are missing in 
|KMSWj . This has no consequence on the validity of the results in this paper, that also 
hold true for the present Hamiltonian H. 

Next we turn to the scattering situation. For simplicity, we restrict ourselves to the 
diatomic case (i.e. M = 2). We still look at the Hamiltionian Pmoi but we want now to 
describe the collision of two ions (or two atoms, or an atom and an ion). It is useful to 
choose a change of variable, that allows a easy description of the system at the beginning 
of the collision process (and another one, to describe the system after the collision). To 
this end, we introduce a cluster decomposition c = {c[,C2} with c'j = {j} Ucj, for j = 1,2, 
and Ci, C2 form a partition of the set {3, ■ ■ ■ , + 2}. At the beginning of the scattering 
process, the particles are gathered in two clusters described by c[ and c'2. Each cluster 
constains a nucleus illustrating the fact that we consider a collision of two ions (and not 
a collision of some electrons with a molecule). Since the motion of the centre of the 
system is not relevant for scattering, we shall remove it. In order to do so, we use the 
particular change of variables in |KMWll IKMW2j . which also allows a good description 
of the scattering processes associated to the decomposition c. 

For k e {1; 2}, denote by \ck\ the number of electron in the cluster k. Its mass is then 
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m'l^ :— rrik + |cfe| and its mass centre is located at: 

Rk ^ (nikZk + ^ . (2.9) 



In particular, the total mass m of the molecule is m = mi + m2 + A?" = m'^ + The 
new variables are, for k e {1; 2}, 



-R := ("^1^1 + rn2Z-2 + ^ 2;^-) , x := i?i - i?2 , := - Zk , for j G . (2.10) 
We set, for y e M=^^, 

<fe)-]^E^.-]4Ew- (2-11) 

As above, this change of variables Cc induces a unitary map Uc by 

{UJ){R-x,y^,--- ,yr,) = Cj{C;\R;x,y,,---,yN)). (2.12) 
Looking at P^oi in the new variables amounts to considering 

2m 

where H' only acts on the variables {x,yi, - ■ ■ , i/tv)- Taking away the motion of the centre 
of mass of the full system again, we keep our attention on H', which is given by 

^HE = -Ei|-E'V„f (2.14) 

k=i jeck 

Q,{x) := Q' + h(x), (2.15) 



^'^""^ " \yj + X - i{y) I " \y, - X + £(y) I + - £(y) I ^^-^^^ 

1 1 

^ .ectTec. \yi-yj + ^- ^(y) I .ectTec 1?/^ - - ^ + ^(?/) I ' 

//^ := H' - I,{x) . (2.18) 

Q'^ stands for the Hamiltonian of separated (noninteracting) clusters while Ic contains all 
extracluster interactions. The electronic Hamiltonian for the nuclear position x is Qc{x). 
The term T^^ is the Hughes-Eckart term in this situation. As usual in scattering theory, 
we should mention the reference (or free) dynamics to which the full dynamics has to be 
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compared (for time tending to — oo). The free dynamics is generated by the Hamiltonian 
H'^ in (12.181) . that is, up to the Hughes-Eckart term T^^, the Hamiltonian of freely moving 
clusters ci and C2. If the scattering process under consideration produces at the end (i.e. 
for time tending to +oo) a cluster decomposition d of the system then the free dynamics 
for large positive times is given by H'^. 

Using Kato's pertubation argument and Hardy's inequality (cf. [Kl IRS2] ). one can show 
that H can be realized as a self-adjoint operator. This means that H can be defined on an 
appropriate subspace of the functions of the variables (xi, ■ ■ ■ ,xm-i,1/i, ■ ■ ■ jUn) and 
the resulting operator is self-adjoint. In fact, this subspace is the domain of the Laplace 
operator in all variables (xi, ■ ■ ■ , xm-i, Hi, - ■ ■ , Un)- The same holds true for Hq. In the 
same way, Pmoi (respectively H') is self-adjoint on the domain of the Laplace operator 
of all variables zi, ■ ■ ■ , zm+n (respectively x,yi, - ■ ■ , i/n)- For fixed x, if we view Q{x) 
(respectively Qc{x)) as an operator on the variables yi, - ■ ■ ,yN, it also has a self-adjoint 
realization by the same argument. If we view Q (respectively Qc) as an operator on the 
variables x,yi, - ■ ■ , yjy, it can be written as a direct integral of self-adjoint operators (cf. 
|SW] ) and therefore has also a self-adjoint realization (cf. p. 284 in |RS4] ) . 

3 Mathematical approach to the Born-Oppenheimer 
approximation. 

Here we want to present the main ideas behind the mathematical treatment of the Born- 
Oppenheimer approximation, which was initiated in [Cl ICDSl ICS] . Since its validity 
should rest on the fact that the nuclei are much heavier than the electrons, one introduces 
a small, positive parameter h related to the nucleon/electron mass ratio. For instance, in 
|KMSWj . the nuclear masses are given by ruk = h'^Xk, where the are of order 1, and, 
in |KMW1[IKMW2] . one uses h'^ = + M^^ (with the notation of Section [2]). Anyhow 
the main point is that h is always sent to 0. This means that the results proved hold 
true for "small enough" h and, most of the time, one has no concrete idea of how small 
h should be. This restriction is of course a drawback for physical or chemical perposes 
but it is useful to understand the small h limit and it gives often correct results when 
compared with the observed behaviour of the physical system. 

Now we come to the main features that ensure the validity of the Born-Oppenheimer 
approximation. Let us consider a normalized, bound state if of energy E of the operator 
H in (12. 5p . We note that = h'^fi'j, where the fi'j do not depend on h. Let us fix the 
nuclear variable x. Typically the spectrum of Q{x) starts with some isolated eigenvalues 
Ai(x), ■ ■ ■ , Aj(x) with finite multiphcity and has above a continuous part ac{Q{x)) (see 
fig. 1 in the diatomic case). Since Q{x) is a self-adjoint Schrodinger operator, we can 
decompose ip{x] ■) in a basis of electronic "eigenvectors" of Q{x): 




(3.1) 
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where the ipj{x; ■) form a basis of true eigenvectors of Q{x) associated to Xj{x) respectively 
and the Ex{x) are spectral projectors of Q{x) with energy A (here (■, ■)y denotes the usual 
scalar product in the y variables). We assume that E belongs to a small open energy 
interval {E^; E+) below the infimum over all x of crc{Q{x)) (as in fig. 1). Since the 
coefficients Ckj in (12.71) contain h"^, the Hughes-Eckart term The is small compared to the 
electronic Hamiltonian Q in (12.81) . To compute E and (f, the idea is that only the part 
of the spectrum of Q{x) less or equal to E+ should be relevant, since the nuclear kinetic 
energy is nonnegative. Define J+ as the largest j < J such that there exists some x with 
\j{x) < In (13.11) . we expect that 



This is the mathematical reading of |BHj . Let n(a;) be the orthogonal projection on the 
first J-)_ levels of Q{x) that is, for an electronic wavefunction ip, 



To implement this idea, it is natural to try as an effective, self-adjoint Hamiltonian an 
operator of the form IIGn where G may be H or Hq. It acts on a total wavefunction 
ip{x,y) as follows: for each x, one projects ip{x, ■) as in (13. 2p then one lets G act on the 
result and, finally, one projects again according to (13. 2p . Thus, the spectral subspaces 
of Q{x) corresponding to energies above E+ are removed. Such an operator is usually 
called an adiabatic operator. We expect that among the eigenvalues in {E^; E^) and 
corresponding eigenvectors of IIGII, there is a good approximation of E and (p. Let us 
try to justify intuitively this claim. 

Assume that {E_] E^) is a small interval around the infimum of the function Ai that is 
attained in some region F (in the diatomic case, Ai only depends on the norm of x and F 
can be a sphere) and that J+ = 1 (for instance (i?_;i?+) = {E^; E^) in fig. 1). We set 
n-'-(x) = 1 — n(x). Using the positivity of the nuclear kinetic energy, H > Q + The thus 



Since The is small compared to Q, one can show via Kato's pertubation theory that 



(p{x,-) = ''^^{ip{x; ■),iljj{x; ■))y4'j{x; ■) + small term . 



J+ 




(3.2) 




{ip{x,-),TYiE<-p{x,-))ydx = 0(/i^) and 



(3.3) 
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where || ■ || is the L^-norm in the variables {x, j/i, ■ ■ ■ , j/tv)- Therefore 

E > j \i{x)\{^{x, ■),iJi{x; ■))y\'^ dx 

+ E+J \\n^ix)ip{x,-)\\ldx + Oih^) , (3.4) 
where || ■ \\y is the L^-norm in the y variables. Since E < E^ and 

1 = = \{v{x,-),i^i{xr))y\'^ dx + / \\U^{x)ip{x,-)\\ldx , 



2 

the integral in (13. 4p and the term \ {^{x, ■),ipi{x; ■))y\ for x far from T should be small. 
In particular, in (13. ip . we should have 

>f{x,-) = x(x)(v9(x; ■),?/'i(x; ■))j^Vi(x; ■) + small term , 

where x is the characteristic function of a neighbourhood of F. The relevant part of ip 
is then its projection onto the electronic level ^i(a;; ■). Since the lower bound in (13. 4p is 
close to E and does not contain the nuclear kinetic energy, the latter must be small. 
Note that this is not sufficient if the interval (-E-; E^) is placed much higher (as in fig. 1). 
Let us now produce a better argument in this more general situation. Using the second 
estimate in (13. 3p . one can show that, close to E and ip respectively, there are an energy 
Eq and a normahzed L^-function ipo such that -ffoV^o = Eoipo. We assume £"0 € (E^; E+). 
In particular, U{Ho-Eo)ipo = and U^{Ho-Eo)ipo = 0. Using U + U^ = 1 and U^ = U, 

iUHoU - Eo)nipo = -UHoU^ipo, (3.5) 

(n^ffon^ - ^o)nVo = -n^Hoiiipo. (3.6) 

Since n(x)-'"n(x) = and n(x) commutes with Q{x), 

M-l M-l 

UHoU^ = -n(5:— A.,)n^ = /.^J:— n[A.^,n]. 

Now /i^ times the commutator [A^.^,,!!] equals 2h{W Xj^){x) ■ hV^^ + h'^{^xj^){x)- Since 

1 2 Z"*^"^ 1 2 

^A^i J til 

and the nuclear kinetic energy remains bounded (this is due to the self-adjointness of 
Hq on the domain of the Laplace operator), the right hand sides of (13.50 and (13.61) are 
0{h) in L^-norm. Since II-^HqII^ > E^, the operator II^HqII-^ — Eq is invertible with 
bounded inverse and (13. 6p shows that H'^ifo is 0{h) in L^-norm. In particular, Hipo is 
almost an eigenfunction of UHoU with L^-norm close to 1 so is close to a true normalized 
eigenfunction of UHqII. Thus ipo (and also ip) should be also close to a normalized 
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eigenfunction of UHqU. Notice that, by (I3.5p . the nuclear kinetic energy in the state Uipo 
is close to the quantity 

- j {11(^0 , {Q{x) - Eo)n(po)y dx 

which has no reason to be small in general. Indeed, if Eq is clearly above the infimum of 
Ai and below the infimum of A2, J+ = 1 and this term equals 

- J (Ai(x) - Eo)\{(^o{x, ■),'4'iix; ■))y\'^ dx 

and it is possible to show that the squared function essentially lives in the "well" : 
{x;Xi{x) <Eo}. 

In the above argument, we used the fact that we can differentiate 11 twice with respect 
to the variable x. This is not obvious at all when one looks at the x-dependence in Q{x) 
(see (12. 8p ) which involves the rather irregular function W in (12. 3p . Thanks to a trick due 
to Hunziker in |Huj . one can prove that x ^ Il{x) is smooth away from the set of the 
nuclear collisions (this is actually sufficient for our argument above since one can show, 
by energy arguments, using the repulsive nature of nuclear interaction, that the wave- 
function is concentrated away from these collisions). This trick is also used in |KMSWj 
and is partially responsible for the technical complications in this paper. The idea here 
is to perform a x-dependent change of variables on the y variables in Q{x) that makes 
the x-dependent singularities in the function W{x,y) in (12. 8 p a;- independent. This can be 
done only locally in x i.e. for x close enough to any fixed position xq (see Lemma 2.1 in 
[J5] for details). 

The regularity of 11 allows the contruction of 7+ globally defined, smooth functions x 1— )■ 
i'ji^x, ■) with values in the electronic functions such that, for each x, ipi{x, ■),■■■ , V^j+(a;, ■) 
is a orthogonal basis of the image of the projection n(x). Because of the possible pres- 
ence of eigenvalues crossing (see fig. 1), it is not always possible to choose these functions 
ipj{x, ■) among the eigenvectors of Q{x). 

We have seen that, up to an error of size 0{h), one can reduce the eigenvalue problem for 
H to the one for UHqU. One can compute explicitely Ili^on in terms of the electronic 
wavefunctions ipi{x, ■), ■ ■ ■ , ipj^{x,-) (see |PSTj ) and again remove terms that are also 
0{h). Writing the approximate eigenstate (p of energy E as 

J+ 

If = '^ej{x)iljj{x; ■) , (3.7) 
i=i 

one ends up, in the diatomic case for simplicity, with the uncoupled equations —h'^A^Oj + 
XjOj = E9j, if the ipj{x; ■) are eigenvectors of Q{x). If the latter is not true (it can be the 
case when crossing eigenvalues occur), one has coupled partial differential equations in 
the X variables for the 6j with coefficients depending on the Xj and on the {ip3i'4'k)y Now, 
if we demand an accuracy of 0{h?), then HHqYI still provides a good approximation but 
less terms can be removed. In particular, one has to keep terms containing the so-called 
Berry connection, i.e. factors of the form {ipjj'Vx ipk)y We see that the variation of the 
iljj (or of n) has to be taken into account. We refer to |PSTj for details. 
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If one wants to improve the accuracy to 0{h^) (or better), one needs to include the 
Hughes-Eckart term The in Q{x) and to replace n(a;) by an appropriate projector Ilg{x) 
which is essentially of the form Il{x) + hlli{x) + h'^Il2{x). Then one uses as adiabatic 
operator n^^i/n^ (cf. [MS21IPST] ). 

In the scattering situation mentioned in Section [21 one can choose the total energy E 
in an energy range {E'^,Et.) like in fig. 2 and take Ilif'II as an effective Hamiltonian. 
In this case, it is important to choose the projections Il{x) as spectral projections of the 
operators Qc{x) + T{j^^ (because the Hughes-Eckart term T^e has no decay in x at infinity). 

4 Review of mathematical results. 

In this section, we present some rigorously proved results on the Born-Oppenheimer ap- 
proximation that illustrate the main ideas developed in Section [31 Since we cannot review 
all results, we selected one from each of the following fields: bound states, resonances, 
scattering process (colhsion) and time evolution. These choices may be detected as abi- 
trary (they reflect the way the author senses the subject) but we try to present results 
with the highest degree of generality. Nevertheless we also comment on other results in 
these fields. At the end of the present section, we add some remarks when symmetries of 
the particles are taking into account. 

Let us begin with the study of bound states of a molecule which was performed in the pa- 
per |KMSWj (previous results were obtained in |Ha31 [Hal] ). One studies the eigenvalues 
of the operator H (cf. (12. 5p ) in the framework introduced in Section [31 In particular, the 
adiabatic operator UHU is used first as an effective Hamiltonian but in a slightly differ- 
ent way. The authors use a so-called Grushin problem and pseudodifferential techniques 
to produce a more accurate effective Hamiltonian E{E) (depending on the sought after 
energy E G (i?_;i?+)), which is a pseudodifferential matrix-operator. E{E) essentially 
corresponds to the operator that defines the (a priori) coupled equations on the 9j we 
mentioned at the end of Section [31 Then E is an eigenvalue of H (essentially) if and 
only if is an eigenvalue of E{E) (cf. Theorem 2.1). Here we mean that, if is a true 
eigenvalue of H, then is an eigenvalue of F{E') where \E — E'\ = 0{h^), for all integer 
A^, and also that, if is an eigenvalue of F{E') then H has an eigenvalue E such that 
\E — E'\ = 0{h^), for all integer N. An explicit but rather complicated, infinite contruc- 
tion produces the operators F{E). For practical purpose, one follows only an appropriate 
finite number of steps of this construction to get an operator Fp{E) such that the above 
errors are 0{W). A concrete example is given as Proposition 1.5, where eigenvalues of H 
in some particular energy range are computed up to 0(/i^/^). 

For diatomic molecules, the authors consider an energy range close to the infimum of 
Ai (like (£'°;S°) in fig. 1). Recall that, for all nuclear positions x, \\[x) is the lowest 
eigenvalue of the electronic Hamiltonian Q(x), which is simple. Actually, one does not 
need to consider the lowest eigenvalue but it is important that it is simple and that the 
rank of the projection H(a;) is always 1. In the mentioned energy range, the eigenvalues 
of H and the corresponding eigenvectors are computed by an asymptotic expansion in 
power of h of WKB type. 
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For polyatomic molecules, the same situation is studied but two cases occur. Recall that F 
denotes the set of nuclear positions x = {xi, ■ ■ ■ , xm-i) ^ ]R^(*'^^^) where the infimum of Ai 
is attained. It is assumed that F is the set of all points (Ox?, ■ ■ ■ , Ox°^_^) where O ranges 
in the set of all orthogonal linear transformations in M.^ and Xq = [x^, ■ ■ ■ ,x%_i) G F. 
One can check if the points x?, ■ ■ ■ ,x^j_^ G lie on a line, or on a plane, or generates 
the whole space M^. The molecule is "linear", "planar", and "non-planar" respectively. For 
a linear or planar molecule, it is shown that, in an appropriate neighbourhood of Ai's 
infimum, there is exactly one eigenvalue of H which is given by a complete asymptotic 
expansion in h. A corresponding eigenvector can also be obtained by such an asymptotic 
expansion. The distance from this eigenvalue to the rest of the spectrum of H is of order 
h^/^. In the non-planar case, two different simple eigenvalues of H are present in the men- 
tioned neighbourhood. The splitting (that is the distance between these two eigenvalues) 
is exponentially small in h. The eigenvalues and the corresponding eigenvectors are again 
given by an asymptotic expansion in power of h. These eigenvectors can be related to one 
another with the help of the reflection ip{x,y) i— )■ (^(— x, —y). 

In the above framework, we mention a modification of the Born-Oppenheimer approxima- 
tion performed in |HJ6t IHJ7t IHJ8] in order to make apparent chemical hydrogen bonds 
in molecules. The main idea is to take the hydrogen mass of order /i^/^, while the mass 
of the heavier atom and the electronic mass are still of order h"^ and = 1, respectively. 
In this setting, one can reduce the eigenvalue problem to an effective one in a similar way 
as in Section 131 However, the authors use a multiscale analysis as in |Ha3| iHai] . 
Next we describe the paper |MM] on the resonances of the operator H in the diatomic case. 
As we shall see, resonances are complex eigenvalues of an appropriate distortion of H. It 
is believed that they give information for the long time evolution of the molecules (scat- 
tering). This link has been done in other (simpler) situations but, to our best knowledge, 
not in the present framework, i.e. for molecules in the Born-Oppenheimer approximation. 
Actually the authors consider the diatomic version of the operator Hq (cf. ( 12. 6p ). that 
is the relevant molecular Hamiltonian without the Hughes-Eckart term. To explain the 
announced distorsion, we need some notation. Let a; : — ?■ a smooth function, which 
is near and equals the identity map (i.e. u{x) = x) for |x| large. For real numbers /i, 
one introduces the transformation defined on total wavefunctions (p{x,y) by 

{U^ip){x,y) = \Jf,{x,y)\^^'^ (fix + nu{x),yi + nu{yi),- ■ ■ ,yN + fiuiyN)) , 

where the function is the Jacobian of the change of variables (x, ?/)(-> (x + ^u){x),yi + 
fiLo{yi), ■ ■ ■ ,yN + l^'^iyN))- It turns out that one can extend to small enough complex 
values of /z. The distored Hamiltonian is given by = U^HqU~^ . Since is unitary 
for real /x, if^ has the same spectrum as Hq on the real line but, for complex /x, the 
continuous part of the spectrum of iJ^ is obtained from the one of Hq by some rotations. 
Furthermore, between the continuous spectra of Hq and if^, eigenvalues of if^ of finite 
multiplicity appear. They actually do not depend on /i and are called the resonances 
of Hq. They are close to the continuous spectrum of Hq., which is responsible for the 
scattering processes governed by Hq. One wants to compute these resonances. To this 
end, one faces an eigenvalue problem as above in [KMSW] but now for the operator H^. 
In [MM] , it is shown that one can adapt the arguments of [KMSW] to show that E is a. 
resonance of a modified version of Hq if and only if is an eigenvalue of an i?-dependent 
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pseudodifferential matrix-operator. Due to a technical difficulty, the authors have to 
smooth out the repulsive, nuclear interaction, leading to the modified version of Hq. Thus 
it is still open whether a similar result holds true for resonances of Hq. Anyhow, the 
imaginary part of these resonances is expected to be of order e"'^/'* in h with c > (as 
shown in jMa2l [Ma8] in a simpler framework). Since the inverse of the imaginary part of 
a resonance is interpreted as the lifetime of the corresponding resonant state, the resonant 
states in the present situation should live on a time scale of order e^^^. 

Next we are interested in the non-resonant scattering (or collision) theory for diatomic 
molecules. We use the framework presented in Section [3] for the operator H' (cf. fl2.13p ) 
but we impose a stronger decay on the pair interactions. We replace at spatial infinity the 
Coulomb interaction | ■ |~^, which has long range, by a short range potential V (essentially 
of the type | ■ 1"^"*^, for some e > 0). To present the result in [KMW2j . we need a short 
review of the short range scattering theory (see details in Section XI.5 p. 75 in |RS3] ) . 
Recall that the free dynamics for large negative time is generated by H'^ in fl2.18p . which 
can be written as 

2 \m\ 1712/ 

We choose an eigenvalue E'' and a corresponding eigenvector ipc of Q'^ + T^p. and consider 
the scattering process that begins for t — )■ — oo by the free motion with kinetic energy 

I \ m\ 7122 ' 

of two clusters that are in the state ipdv)- The initial state is described by a wavefunction 
T{x)ipc{y), where r is a "nuclear" wavefunction and ipdv) is actually the product of an 
electronic wavefunction of the cluster ci by one of the cluster C2- If we forget about 
the Hughes-Eckart contribution T^g, these latter wavefunctions represent respectively an 
electronic bound state in the atom/ion ci and another in the atom/ion C2- It turns out 
that we can find a wavefunction (p+{x,y) such that, for t — )■ — oo, the real evolution of 
ip+{x,y) is close to the free evolution of T{x)ipc{y), that is 

||e-^*^V+ - e-'''''^T{x)My)\\ ^ 0. 

We denote by fi+ the operator T{x)ipc{y) ^ v^+. Similarly, the final state is descbribed 
by t' {x')ipd{y') corresponding to some cluster decomposition d (with a priori different 
coordinates) and one can find a wavefunction (f_{x',y') such that, for t — i- +oo, 

||e-'*^V- - e-'*^^r'(x')V^rf(l/')|| ^ 0. 

We define Q^{t' {x')ip'^{y')) = Lp_ and we can check that t' {x')ip'^{y') can be recovered from 
(p-{x',y') by application of the adjoint Q*_ of Thus the operator S = sends 
the initial state to the final state of the scattering process. It is the scattering operator 
for this process while f2+ and Q- are the wave operators. The strange sign convention for 
the wave operators can be interpreted in the following way: e'^*^^ Qj^T{x)tpc{y) represents 
the future (+) evolution for the interactive dynamics (defined by H') of the free state 
T{x)ipc{y)- When c = d, t = t' , and ipc = we have an elastic scattering process. When 
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c = d but ipc and ^/'^ are orthogonal, the inelastic process corresponds to a change of 
electronic level in the cluster (an excitation of an electron in ci for instance). When c ^ d 
but ci and di contain the same nucleus and so do C2 and d2, an electron at least has moved 
from one nucleus to the other. We can also consider the case where c is as above while di 
contains the two nuclei and d2 only electrons (for instance, two ions form a molecule and 
loose some electrons). Among the inelastic processes we just described, the two last ones 
might be interesting for Chemistry. The above construction of a scattering operator can 
be done for all possible cluster decompositions c and d and the collection of S operators 
completely describes the possible scattering processes. The same construction can also 
be performed for molecules with more than 2 nuclei with a richer family of processes of 
chemical interest. 

We point out that a scattering theory exists for long range interaction (like the Coulomb 
one). Essentially, one has to modify the construction of the wave operators, which become 
technically more involved. 

We come back to the situation studied in |KMW2] . that is for a diatomic molecule with 
d = c given as in Section [2] but with short range interactions. We choose an energy range 
{E'i,E'^) as in fig. 2. In particular, it is above the infimum of the spectrum of H'^ in 
(12.181) thus, by the HVZ Theorem (see Theorem XIII. 17 p. 121 in |RS4j ) . this energy 
range is included in the continuous part of H' but might contain eigenvalues. We focus 
on scattering processes with total energy E G (-E1, E1_). In view of Section [3l we replace 
E+ by i?^ and consider eigenvalues \j{x,h) of Qc{x) + T^p^ that are somewhere less or 
equal to E'^. One then constructs the associated projections Il{x,h) and consider the 
adiabatic operator UH'U. Let fl"^ be wave operators associated to the decomposition c 
as above such that the electronic energy E'^ of the initial state satisfies E'^ < E'L [E^ is 
close to some eigenvalue E'^{{)) of Q'^, see fig. 2). Now one can also compare the dynamics 
of YlH'Ii and the free dynamics generated by H'^ and construct so-called adiabatic wave 
operators The goal is to show that 

- ^]^°|| = o{h) , (4.1) 

when the wave operators act on wavefunctions with energy in {E'L, E'^). To this end, one 
needs an important assumption, the non-trapping condition, on the classical mechanics 
generated by the classical Hamilton functions hj{q,p) = + Xj{q,0) (with h = and 
for the above selected eigenvalues A^) at energy in {Et, E'^). This non-trapping condition 
says that all classical trajectories of energy E G {E'i,E'^) for any Hamilton function hj 
go to spatial infinity in both time directions. It implies the absence of eigenvalues in 
{Et,E'^) and prevents resonance phenomena. In fig. 2, this assumption is satisfied. 
Under the additional assumption that only the simple eigenvalue Ai is somewhere less or 
equal to E^, (in particular the image of H(a;, h) is always of dimension 1), the approxima- 
tion (14. ip is proved in |KMW2] . An important step in the proof is to etablish an appro- 
priate estimate on the resolvents C\M9zi— )'(2; — H')~^ and C\]R92;i— )-(z — UH'lVj'^ 
of H' and UH'U respectively and this is done for long range interactions (in particular 
for the Coulomb one). Because of the additional assumption, only elastic scattering is 
covered. In this framework, we mention the papers [J2] on the scattering operator and 
|JKW] on scattering cross-sections. 

If one removes the above additional assumption, one can obtain the approximation (14. ip 
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but under the condition that the eigenvalues Xj do not cross (see [JTj). In this situation, 
a similar approximation holds true for scattering cross-sections and it can be shown that 
the inelastic scattering is disadvantaged compared to the elastic one (see [J3]). In the 
simplified framework of a Schrodinger operator with matrix potential, it is even shown in 
|BM] that the inelastic scattering is exponentially small in h. Therefore, to study it, we 
have to accept eigenvalues crossings and we need to control their effect on the scattering. 
As mentioned before, the projection 11 is still smooth but the eigenvalues Xj might be only 
continuous and the corresponding eigenvectors ipj might be discontinuous at the crossing. 
In this situation, we mention the work by [FR] on Schrodinger operators with matrix po- 
tential and for a special case of crossing (crossing at just one point), where the resolvent 
estimates mentioned above are derived. For some types of crossing, the Xj and ipj are 
smooth and one can prove the same result (see p4| IDFJj ). This is the case for diatomic 
molecule thanks to the radial symmetry of the molecule with respect to the x variable. In 
the work in progress [JS], one uses this property to get the resolvents estimates and also 
the approximation (14. ip for diatomic molecules. 

Now we come to the last field we wanted to consider, namely the time evolution of 
molecules in the Born-Oppenheimer approximation, and present results obtained in |MS2] . 
We consider again the operator H in (12. 5p but we look for an approximation of the evolu- 
tion operator e"**^/'' (i.e. the molecular evolution on a time scale 1/h). As in Section [3l 
the authors choose a certain energy range (like {E_; E^) in fig. 1) and construct a better 
projection II^, starting from the operator 11 adapted to this energy range. The estimate of 
the commutator [H, U] = 0{h) is improved in this way in the estimate of [H, Ug] = 0{h^), 
for all integer p. With the help of Ug, the authors introduce a map W that transforms 
wavefunctions ip{x,y) for the full molecule into wavefunctions in x only but with values 
in the L-dimensional vectors {L being the constant dimension of the image of the Il{x)). 
This map replaces the electronic wavesfunctions, that live in an infinite dimensional space, 
by a finite number of degrees of freedom, namely the coordinates of the L-dimensional 
vectors. We point out here that no restriction on the number of nuclei is required and 
that eigenvalue crossings are allowed. There exists a L x L-matrix operator A acting on 
the range of W such that, for a large class of initial states (f, for all integer p, the time 
evolution of cp is given by 

^-^tH/h^ = W*e-''^^^Wip + 0{{l + \t\)hP) , (4.2) 

where t ranges in some bounded, p- and /i-independent interval. So, to compute a good 
approximation of the time evolution of ip, one first lets W act, then follow the evolution 
of yV(p generated by A (a much simpler evolution) and then lets the adjoint of W act. 
The operator A is obtained by an infinite but explicit procedure. If one accepts an error 
of size 0{{1 + \ t\)hP), for a fixed p, one can replace A by an operator Ap which is obtained 
by a finite procedure. 

As a consequence of the previous approximation, the authors derive for L = 1 a rather 
precise description of the time evolution of coherent states (which are probably the sim- 
plest states), completing in this way previous results of this type (for instance |HaltlHa5] ). 
The assumption L = 1 prevents eigenvalue crossings. For the time evolution of coher- 
ent states, the effect of eigenvalue crossings was studied in |Ha7j . Even for these states. 
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this effect is complicated in general and another approach was followed by considering 
so-called avoided crossings (see |HJH [HJ2] ). Instead of having a crossing of the electronic 
eigenvalues Ai and A2, one assumes that, for some particular nuclear position, the nonzero 
difference Ai — A2 is small (with an appropriate size compared to h). This approach avoids 
the technical difficulties carried by true crossings but allows inelastic phenomena (like the 
transfer of a wave packet from the electronic level Ai to A2). In a simplified framework 
(compared to the molecular setting) but for the time evolution through true eigenvalue 
crossings, we mention |FG] in a special case where the Xj and the eigenvectors ipj are not 
smooth and jPFJj where the latter are smooth. In |FG] a Landau-Zener formula plays 
an important role. In |DFJj . although the coupling of the smooth crossing eigenvalues 
vanishes formally at = 0, a coupling effect between them is proved in a very special 
situation, that should be unphysical. Finally we quote the paper [TWj where the Born- 
Oppenheimer approximation for the time evolution of molecules coupled to a quantized 
radiation field is analysed. 

We end this section with some comment on the symmetries of particles. First one should 
take into account that the electrons are fermions and consider only antisymmetric elec- 
tronic wavef unctions. Second, if the molecule contains two identical nuclei for instance, 
one should restrict the nuclear wavefunctions to the ones that are symmetric with respect 
to the exchange of these two nuclei. In principle, such constraints can be included in a 
mathematical framework but, in practice, this has not been done. Let us give some expla- 
nation for this. Including these symmetries amounts to let act the operators on smaller 
Hilbert spaces. So if one can perform the approximation in the full Hilbert space, it is also 
valid on a smaller one. However, the electronic symmetry could change the spectrum of 
the electronic Hamiltonian (the eigenvalue A3 could be absent or its multiplicity could be 
lowered) but this would change essentially the imput of the above mathematical treatment 
and not the core of the approximation. Taking into account the nuclear symmetry could 
give finer results but this would be hidden in the properties of the adiabatic operator or 
the effective operator derived by the mathematical Born-Oppenheimer approximation. Up 
to now, it seems that there was no clear motivation from the mathematical point of view 
to include symmetries; thus it was natural to avoid them and the technical complications 
they carry. 

5 Conclusion. 

We have presented the essential structure of the mathematical justification of the Born- 
Oppenheimer approximation and tried to illustrate it on concrete results on bound states, 
on the time evolution, and in scattering theory. In particular, we have seen that the 
basic idea is close to the point of view presented in |BHl ISWj and consists in writing the 
full Hamiltonian as the sum of the nuclear kinetic energy, of an electronic Hamiltonian, 
and of comparatively smaller terms, mimicking in this way the usual framework for the 
well-developed semiclassical analysis. Indeed, taking the favorite example of this analysis, 
namely the semiclassical Schrodinger operator —h'^A^ + V{x), the nuclear kinetic energy 
stands for the semiclassical Laplace operator — /i^A^, while the electronic Hamiltonian 
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plays the role of the potential V. We have explained how the full Hamiltonian can be 
approximated by a so-called adiabatic operator, the construction of which essentially rests 
upon the electronic Hamiltonian (or clamped-nuclei Hamiltonian). Even the construction 
of the refined projection H^, which leads to a very accurate approximation, completely 
depends on this Hamiltonian. We point out that our intuitive argument to compute an 
eigenvalue and an eigenvector of the full Hamiltonian (the operator H), up to an error 
0{h), actually leads to Born-Huang's proposition of approximated eigenvalue and eigen- 
vector (see fl3.7p ). Born-Huang's approach is legitimate but not very accurate. To go 
beyond, as we mentioned, one needs to take into account the variation of the electronic 
Hamiltonian with respect to the nuclear variables. When we look for an eigenvalue close 
to the groundstate energy (which is close to the infimum of the lowest electronic eigenvalue 
Ai), we have seen that the nuclear kinetic energy is small, as a consequence of this closeness 
and not of the large size of the nuclear masses. In particular, the original computation in 
|B0] is legitimate. The situation is different for higher energy but it can be handled with 
the help of semiclassical analysis (see |KMSW] ) , as explained in Section HI Concerning 
the scattering (or collision) theory and the time evolution of molecules, we reviewed some 
results and pointed out the main difficulty, namely the control of eigenvalue crossings. In 
particular, this difficulty hinders the treatment of chemically relevant situations but we 
stressed that some progress was made. Letting h tend to instead of keeping its physical 
value is essential in all the above mathematical works but, as we noticed, it might be 
inappropriate in some physical or chemical situations. 

In |SWj . the authors subscribed to Lowdin's impression (expressed in [L]), that it might 
be difficult to extract from the molecular Hamiltonian the concrete realization of chem- 
ical concepts like isomerism, conformation, chirality. Probably, they are right but the 
situation is perhaps not hopeless. We pointed out the papers |HJ6t IHJ71 IHJ8] that try 
to describe hydrogen bonds. In the paper [JKW] . it was proved that some symmetries 
in the ion-atom scattering influence the leading term of scattering cross-sections in the 
Born-Oppenheimer approximation. The techniques used in [KMSWj tells us that, near 
the minimum of a nondegenerate electronic eigenvalue Xj, one can flnd a bound state of 
the molecule with low nuclear kinetic energy. In this state, the nuclei vibrate near an 
equilibrium position, located where the minimum is attained. If one can compute (nu- 
merically) this position, one gets the nuclear structure of this bound state (internuclear 
distances, symmetries). Because of computational error, it might be difficult to check if 
the molecule is planar or not. By light excitation, one can measure the difference between 
the molecular energies, that are the two closest levels to the minimum of A-,. If the differ- 
ence is "very" small, then the molecule in this state is not planar and if the difference is 
big "enough", then it is planar, thanks to [KMSWj . Of course, these examples are limited 
from the chemical point of view but show that simple properties of molecular structure 
can be extracted from the molecular Hamiltonian. We also emphasize that there exist 
tools, like the theory of (co-)representations, to take into account symmetries of molecules. 
An example of such use in the molecular context is provided in |Ha7] . 
We emphasise that, in the mathematical treatment of the Born-Oppenheimer approxima- 
tion, the nuclei are always considered as quantum particles. The use of clamped nuclei 
is just a tool to construct an appropriate effective Hamiltonian but the latter is a quan- 
tum, nuclear Hamiltonian with restricted electronic degrees of freedom. It seems that the 
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Born-Oppenheimer approximation in Chemistry often reduces to the computation of the 
molecular potential energy surfaces (that are essentially the energy surfaces corresponding 
to our classical Hamiltonians hj) and to classical motion of the nuclei on these surfaces. 
Except for the special situation studied in |BQ] (concerning the ground state), a classical 
treatment of the nuclei is not justified, as already pointed out in |SW] . This does not 
mean that classical behaviours of the nuclei do not emerge. On the contrary, it is a gen- 
eral fact that semiclassical situations (like the molecular one with small h) are strongly 
connected to classical features and often produce effects that ressemble classical ones at 
the macroscopic level. For instance, we have seen that the classical Hamiltonians hj play 
a role in the scattering situation and explained above that, under some circumstances, 
the nuclei can be viewed as classical particles vibrating near an equilibrium position. The 
method of potential energy surfaces can be applied to find some molecular excited energy 
near the bottom of an electronic potential well but may not capture all of them. So we 
subscribe to the warning addressed to Chemists in [SWj . 

To the presentation of the Born-Oppenheimer approximation in [SWj . we essentially added 
the notion of adiabatic Hamiltonian as a central tool and insisted on the semiclassical 
structure of the approximation that is close to the model —h^/S.^ + V{x). Our optimistic 
view of future mathematical developments actually relies on the power and the diversity 
of techniques provided by the semiclassical analysis. 

The actual mathematical treatment of the Born-Oppenheimer approximation for molec- 
ular systems is expressed in a rather involved language and provides a theoretical infor- 
mation on such systems, that might be considered as unsatisfactory from the physical or 
chemical point of view. We tried to make it accessible to a large audience and to show 
that, despite the real difficulties it has to face, it could be improved, taking more and 
more into account physical and chemical preoccupations. 
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